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Abstract 

The  dynamical  aspects  of  solid-solid  phase  transformations  are 
studied  within  the  framework  of  the  theory  of  thermoelasticity.  A 
problem  of  the  Riemann  type  for  a  one-dimensional  bar  undergoing 
an  adiabatic  process  is  analyzed.  It  is  shown  that  by  imposing  a  ki¬ 
netic  relation  and  a  nucleation  criterion  it  is  possible  to  single  out  a 
unique  solution.  This  extends  to  the  thermomechanical  case  results 
previously  found  in  a  purely  mechanical  context. 
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1  Introduction. 


Solid-solid  phase  transitions  can  occur  in  materials  by  a  change  in  stress 
or  temperature.  At  a  macroscopic  level,  within  the  framework  of  thermo¬ 
elasticity,  materials  that  exist  in  different  solid  phases  can  be  modeled  by 
a  nonconvex  Helmholtz  free  energy  density  that  depends  smoothly  on  the 
deformation  gradient  and  the  temperature.  At  some  fixed  temperature,  the 
Helmholtz  potential  exhibits  multiple  wells  as  a  function  of  the  deformation 
gradient.  Each  well  is  associated  to  a  distinct  solid  phase. 

In  a  one-dimensional  setting,  Abeyaratne  &:  Knowles  [4]  analyzed  an 
initial  value  problem  of  the  Riemann  type  for  a  bar  made  of  a  special  ther¬ 
moelastic  material  that  can  sustain  multiple  phases.  It  was  found  that  the 
solution  to  the  adiabatic  Riemann  problem,  which  involves  initial  data  from 
two  distinct  phases,  is  not  unique.  Because  of  the  nonconvexity  of  the  the 
Helmholtz  potential,  the  structure  of  the  solution  involves  two  types  of  dis¬ 
continuities,  viz.,  shock  waves  and  phase  boundaries  (see,  e.g.,  Abeyaratne 
&  Knowles  [1]-[4],  Shearer  [7],  Slemrod  [8],  Truskinovsky  [10]).  The 
solution  can  be  determined  up  to  the  phase  boundary  velocity. 

Based  on  the  formalism  of  non-equilibrium  thermodynamics  (see,  e.g., 
C ALLEN  [5],  Truesdell  [9])  and  classical  concepts  used  in  materials  science 
(see,  e.g.,  Christian  [6]),  Abeyaratne  &  Knowles  propose  the  existence 
of  an  additional  piece  of  constitutive  information,  a  kinetic  relation,  which 
applies  at  subsonic  phase  boundaries.  The  kinetic  relation  relates  the  phase 
boundary  velocity  to  the  driving  traction,  which  is  associated  to  a  jump  in  en¬ 
tropy  across  the  phase  boundary.  With  this  additional  relation,  it  is  possible 
to  single  out  a  unique  solution  for  the  Riemann  problem  (see  Abeyaratne 
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&  Knowles  [1]-[4]). 

The  purpose  of  the  present  analysis  is  to  investigate,  within  the  context 
of  the  adiabatic  theory,  the  Riemann  problem  involving  initial  data  from  the 
same  phase.  The  results  are  qualitatively  different  from  those  in  [4]  in  the 
sense  that,  in  order  to  recover  uniqueness,  a  nucleation  criterion  has  to  be 
enforced.  This  nucleation  criterion  signals  the  onset  of  a  phase  change  and 
generalizes,  for  the  thermoelastic  case,  the  results  found  in  [3]  in  the  context 
of  a  purely  mechanical  theory. 

In  Section  2,  the  governing  field  equations  and  the  corresponding  jump 
conditions  for  an  adiabatic  process  are  obtained  from  balance  principles.  A 
particular  thermoelastic  material  is  introduced  in  Section  3.  The  field  equa¬ 
tion  and  jump  conditions  are  specialized  for  this  material  in  Section  4.  An 
initial  value  problem  of  the  Riemann  type  for  a  one-dimensional  bar  is  formu¬ 
lated  in  Section  5.  It  is  found  that  there  exist  two  different  kinds  of  solutions. 
One  kind  involves  no  phase  transition,  whereas  the  other  kind  involves  a  two- 
parameter  family  of  solutions  containing  two  propagating  phase  boundaries 
(the  phase  boundary  velocities  acting  as  parameters).  In  the  latter  case,  a 
kinetic  relation  is  used  to  obtain  the  phase  boundary  velocities  and  hence 
the  solution  is  fully  determined.  It  is  also  found  that  for  some  range  of  the 
initial  data,  it  is  possible  to  have  either  kind  of  solution  (i.e.,  with  or  without 
phase  boundaries) .  Therefore,  a  nucleation  criterion  is  enforced  to  single  out 
the  unique  solution  to  the  problem  from  the  two  different  kinds.  Finally,  in 
Section  7,  the  connection  between  the  Riemann  problem  with  uniform  data 
and  the  nucleation  criterion  is  investigated. 
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2  Prelimin£tries. 


Consider  longitudinal  motions  of  a  bar  with  uniform  cross  section  that  oc¬ 
cupies  the  interval  [— L,  L]  in  an  unstressed  reference  configuration  during 
a  time  interval  [tot^i]-  Assume  that  the  referential  mass  density  p  is  uni¬ 
form.  Let  (x,t)  6  [— L,  L]  X  [to,ti]  be  a  point  in  the  reference  configura¬ 
tion  that  is  mapped  to  the  position  y{x,t)  =  x  +  u{x,t),  where  u  is  the 
displacement.  It  is  assumed  that  u  €  C'°([— L,  L]  x  [to,^i])  and  piecewise 
C^{[—L,  L]  X  [to,  ti]).  Let  7  =  and  v  =  Uthe  the  strain  and  particle  veloc¬ 
ity  respectively  (the  subscript  refers  to  partial  differentiation).  The  restric¬ 
tion  —1  <  7(x,  t)  ,Vx,  t,  guarantees  that  the  deformation  y{',t)  is  invertible 
at  each  t.  Let  a  be  the  stress,  e  the  internal  energy  per  unit  referential  mass 
and  77  the  entropy  per  unit  referential  mass.  Assume  that  a,  e  and  77  are 
piecewise  C'^([— L,  L]  x  [to,^i])-  The  balance  of  linear  momentum,  balance 
of  energy  and  the  Clausius- Duhem  inequality  (dissipation  inequality)  for  an 
adiabatic  thermomechanical  process  in  a  one-dimensional  bar  are 


m::  = 

d 

rX2 

(1) 

dt . 

L 

pvdx  , 

= 

d 

rx2 

/  1  \ 

(2) 

dt . 

J  Xl 

pU+2^  Ida:  , 

d 

pX2 

(3) 

m- 

dt . 

J  X\ 

pridx  >  0  , 

Vt  €  and  V[xi,X2]  Q  [— L, L].  Here,  r{t)  is  the  entropy  production 

rate  in  the  interval  [xi,  X2].  Localization  of  the  balance  laws  (l)-(3)  at  points 
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where  the  fields  are  smooth  provides 


cTx  =  m , 

(4) 

(crr;)x  =  p  , 

(5) 

Vt>0  . 

(6) 

The  compatibility  equation  is 

Vx  =  It  ■ 

(7) 

Equation  (5),  with  the  use  of  (4),  can  be  expressed  as 

CT'Tt  =  P^t  • 

(8) 

Suppose  that  there  is  a  point  x  =  s{t)  in  [xi,X2]  at  which  some  (or  all)  fields 
are  discontinuous.  Localization  at  x  =  s{t)  of  the  global  balances  (l)-(3)  and 
the  compatibility  equation  (7)  provide  the  corresponding  jump  conditions, 
viz., 


M  +  s[7l  =  0  , 

(9) 

[a]  +  ps|u]  =  0  , 

(10) 

+  psfe  +  =  0  , 

(11) 

IvP  <  0  , 

(12) 
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where  s  =  ds/dt  is  the  speed  of  propagation  of  the  discontinuity  and,  for  any 
function  g, 

l5l  =9'^  -9~  ,  9'  =  lim  g{x,  t)  ,  g-^  =  lim  g{x,  t)  . 

Discontinuities  are  classified  into  two  types,  viz.,  classical  shock  waves  and 
phase  boundaries.  Shock  waves  are  related  to  discontinuities  where  the  ma¬ 
terial  on  each  side  is  in  the  same  phase,  whereas  phase  boundaries  refer  to 
discontinuities  where  different  material  phases  exist  on  each  side.  In  the  adi¬ 
abatic  theory,  a  shock  wave  whose  Lagrangian  velocity  s  is  zero  is  referred 
to  as  a  contact  discontinuity  (i.e.,  the  material  particles  on  each  side  of  the 
discontinuity  are  the  same  at  all  times).  Specific  jump  conditions  in  each 
case  axe  given  in  Section  4. 

Let  the  Helmholtz  potential  be  given  by  ^  =  •0(7,^).  This  potential  is 
related  to  the  internal  energy  through 


Ip  =  e  —  Or]  .  (13) 

For  a  classical  thermoelastic  material,  the  stress  a  and  the  entropy  r]  are 
given  by 


a  =  ai'f,  e)  =  fnpiil,  Q)  , 

(14) 

7]  =  6)  =  -'ipei'r,  0)  . 

(15) 

The  isothermal  elastic  modulus  p,  the  specific  heat  at  constant  strain  c  and 
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the  coefficient  of  thermal  expansion  a  are  defined  by 


_  \ 

fj,  =  ^(7,  e)  =  0-^(7,  e)  =  fyipjyij,  e) , 

c  =  c(7,  e)  =  Orieij,  0)  =  -Bxjjeei'y,  0)  , 

>  (16) 

The  isothermal  sound  speed  is  defined,  when  //  >  0,  as 

„  =  ah. 

(17) 

Proper  phase  boundaries  travel  at  subsonic  velocities  (see  also  Truski- 
NOVSKY  [10]).  Using  equations  (13),  (14)  and  (15)  in  equation  (8)  provides 

an  alternative  expression  for  the  energy  equation,  i.e.. 

Vt  =  0  . 

(18) 

It  follows  that  the  dissipation  inequality  (6)  is  trivially  satisfied  at  regular 
points.  The  rate  of  entropy  production  for  a  segment  [xi,  2:2]  of  the  bar  which 
contains  a  propagating  discontinuity  at  x  =  s{t)  can  be  expressed  as 

T{t)  =  rb{t)  +  T,{t) , 

(19) 

where 

pX2 

^b{t)  =  /  PVtdx  , 

J  Xi 

(20) 

^s{t)  =  -pIv]s  . 

(21) 
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Here,  Tb  represents  the  bulk  entropy  production  and  corresponds  to  the 
entropy  production  due  to  the  moving  discontinuity.  Equations  (18)  and  (20) 
imply  that 


Fft  =  0  ,  (22) 

hence  the  entropy  production  for  a  thermoelastic  material  under  an  adiabatic 
process  occurs  solely  because  of  the  presence  of  a  moving  discontinuity.  Based 
on  this  entropy  production,  define  the  driving  traction  as 

f=-PlvW)-  (23) 

For  a  discussion  of  the  notion  of  driving  traction,  see  Abeyaratne  & 
Knowles  [1]  and  Truskinovsky  [10].  The  velocity  jump  can  be  elimi¬ 
nated  from  (11)  by  using  equations  (9)  and  (10),  i.e., 

(pfe]  -  (<7)l7l)s  =  0  ,  (24) 


where,  for  any  function  g, 

{9)  =  \{9^  +  9~)  • 

From  the  energy  jump  condition  (24),  when  i  0,  it  follows  that 

pH  =  (<r)|7l  . 

Moreover,  since  {rjO}  =  (77)[^1  +  M(^),  then  (13),  (23)  and  the  above  relation 
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provide  the  following  equivalent  expression  for  the  driving  traction: 


/  =  pm  -  (^)l7l  +  p{v)P]  ■  (25) 

For  a  thermoelastic  material,  (14),  (15)  and  (25)  give 

/  =  p{m  -  (V'7)l7l  -  m)M}  ■  (26) 

In  terms  of  the  driving  traction,  the  dissipation  inequality  (12)  can  be  ex¬ 
pressed  as 


fs>0.  (27) 

The  driving  traction  plays  an  important  role  in  the  kinetic  relation  and  the 
nucleation  criterion  as  explained  in  Section  6. 

3  Thermoelastic  material. 

Abeyaratne  &  Knowles  [2]  proposed  an  analytical  model  of  a  thermo¬ 
elastic  material  that  can  exist  in  different  solid  phases.  The  model  is  simple 
yet  it  incorporates  the  fundamental  characteristics  of  a  phase-changing  mate¬ 
rial  and  allows  to  obtain  results  in  specific  problems.  A  phase  diagram  of  this 
material  in  the  7,0-plane  is  shown  in  Figure  1.  For  temperatures  below  a  crit¬ 
ical  temperature  6c,  the  material  can  exist  in  either  a  low  strain  phase  Pi  or 
a  high  strain  phase  P3.  These  phases  are  thermodynamically  metastable  and 
are  separated  by  an  unstable  phase  P2.  Above  the  critical  temperature  the 
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material  can  only  exist  in  a  stable  phase  P.  Throughout  this  analysis,  only 
transformations  from  or  to  the  low  and  high  strain  phases  are  considered. 
Moreover,  for  a  fixed  temperature  below  the  transformation  temperature  Ox, 
the  absolute  minimum  of  the  Helmholtz  potential  corresponds  to  the  high 
strain  phase  and  for  a  fixed  temperature  between  the  transformation  tem¬ 
perature  and  the  critical  temperature,  the  absolute  minimum  corresponds 
to  the  low  strain  phase.  A  thorough  description  of  the  thermomechanical 
characteristics  of  this  material  can  be  found  in  [2].  The  boundaries  between 
the  different  phases  are  given  by 


7m(^)  =  7c  +  A/'(0  —  9c)  , 
7m  (^)  =  7c  + 1^(9  -  9c)  , 


(28) 


where  7c  >  0,  9c  >  0,  M  and  m  are  constants.  The  expression  for  the 
Helmholtz  potential  is  given  in  each  phase  by 


-  —y{9  -  9t)  -  c9  log  ( in  Pi, 

2p  p  \^T J 

^(7  _  - ^(7  -  7t)(^  -  9t) 

2p  p 

-c9  log  ( -f  ^{9  -  9t)  in  P3, 

\UT J 


(29) 


where  7^  is  the  transformation  strain  and  Xt  is  the  latent  heat  at  the  trans¬ 
formation  temperature.  The  remaining  parameters  were  defined  in  Section  3 
and  are  assumed  constant.  The  expression  in  the  unstable  phase  P2  is  not  rel¬ 
evant  here.  According  to  the  model  developed  in  [2] ,  the  material  parameters 
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must  satisfy  the  following  restrictions: 


7r  >  (M  —  m)6c  >  0  , 

7c  =  ^  +  +  ’^)(^c  -  ^t)  , 

M  +  m=  +  2a  . 


(30) 


The  first  restriction  guarantees  that  there  is  no  overlap  between  the  two 
metastable  phases  in  the  7,0-plane  (hence,  the  stress  is  uniquely  determined 
for  a  given  temperature  and  strain).  Restrictions  (30)2,3  are  related  to  the 
fact  that  the  metastable  phases  and  the  unstable  phase  coincide  at  the  critical 
point  (7C)^c)-  Prom  (14)  and  (29),  the  stress  response  function  is  given  by 


/i7  -  afj,{9  -  9t) 


<7(7,0)  =  { 


in  Pi, 


y  /i(7  —  7t)  —  a/r(0  —  9t)  in  P3. 


(31) 


Observe  that,  for  a  fixed  temperature,  the  stress-strain  relation  is  linear  in 
eaoh  phase.  Hence,  this  material  will  be  referred  to  as  the  trilinear  material. 
The  entropy  response  function,  from  (15)  and  (29),  is  given  by 


in  Pi, 


Vil,  { 


\  1  /  ^  \  .  T, 

— (7-7r)  +  clog  —  M-C-—  m  P3. 

,  p  \“t J 


(32) 
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It  is  convenient  to  introduce  a  set  of  nondimensional  parameters  for  the 
trilinear  material.  Define  the  following  parameters: 


The  variable  S  might  be  viewed  as  a  “normalized”  strain. 


(33) 


4  Field  equations  and  jump  conditions. 

Henceforth,  for  simplicity,  the  coefficient  of  thermal  expansion  is  taken  as 
zero.  Consequently,  the  mechanical  and  thermal  effects  are  decoupled  in  the 
differential  equations.  Nevertheless,  a  connection  between  the  strain  and  the 
temperature  remains  in  place  via  the  jump  conditions.  For  the  thermoelastic 
material  introduced  in  Section  3,  at  points  where  the  fields  are  smooth,  (4), 
(7)  and  (8)  provide,  with  a  =  0, 

7t-Vx=0  ,  (34) 

Vt  -  (^7x  =  0  ,  (35) 

=  0  .  (36) 

The  stress-strain  relation,  with  a  =  0,  is  shown  in  Figure  2.  The  jump 
conditions  (9)- (12)  at  a  point  of  discontinuity  where  both  sides  are  on  the 
same  phase,  when  a  =  0  and  using  the  nondimensional  quantities  (33),  are 
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given  by 


|(5]v  +  |t;l  =  0  , 

W  +  I^lv  =  0  , 
irjv  =  0 , 

/r+\ 

log  ( ^  J  V  <  0  . 


These  jump  conditions  are  equivalent  to  either 


V 


v[51  +  M 

in 


±1 , 
0, 
0, 


or 


V  =  0, 
<  {Sj  =  0, 
Jt5j  =  0. 


(37) 


(38) 


The  first  case  corresponds  to  a  shock  wave  whereas  the  second  is  a  contact 
discontinuity.  Contact  discontinuities  are  stationary  in  the  reference  config¬ 
uration  (i.e.,  in  the  Lagrangian  sense).  If  the  high  strain  phase  is  on  the  right 
of  the  phase  boundary,  then  the  jump  conditions  (9)-(12),  when  a  =  0  and 
in  nondimensional  form,  are  given  by 


vM  -i-  lu]  =  0  , 

|5]  -  1  -h  vful  =  0  , 
{|ri-((i>-^)-ir}v  =  o, 


(39) 
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If  the  high  strain  phase  is  on  the  left  of  the  phase  boundary,  then  the  jump 
conditions  (9)- (12)  are  given  by 


vM  +  |u]  =  0  , 

I5l  +  l  +  v[t;l  =  0, 

{[T1+((7)-5)+'t|v  =  0, 


(40) 


5  Riemann  problem. 


Consider  the  Riemann  problem  corresponding  to  an  infinitely  long  bar  com¬ 
posed  of  the  thermoelastic  material  described  in  Section  3.  This  problem 
can  be  formulated  as  follows:  Find  functions  5{x,t),  v{x,t)  andT{x,t)  such 
that  equations  (34)-(36)  are  satisfied  at  points  where  the  functions  are  suffi¬ 
ciently  smooth  and  the  corresponding  jump  conditions  (38)- (40)  are  satisfied 
at  points  where  the  functions  are  discontinuous,  subject  to  the  following  initial 
conditions: 


(5(x,0)  ,  v{x,0)  ,  T{x,0) 


1  Tl  for  —00  <  X  <  0  , 
Sr  ,vr  ,Tr  for  0  <  x  <  oo  . 


(41) 


14 


It  is  assumed  that,  initially,  the  bar  is  in  the  low  strain  phase,  in  which  case 


h  G  (— , 

Sr  ^  (— • 

The  structure  of  the  solution,  as  shown  by  Abeyaratne  &  Knowles  [4], 
must  necessarily  involve  either  no  phase  boundaries  (in  which  case  all  par¬ 
ticles  of  the  bar  remain  in  the  original  phase  at  all  times)  or  two  phase 
boundaries  moving  in  opposite  directions  (in  which  case  the  particles  jump 
to  the  high  strain  phase).  In  the  latter  case,  as  shown  below,  there  is  a 
two-parameter  family  of  solutions. 

5.1  Solution  with  no  phase  boundary. 

The  solution  with  no  phase  boundary  involves  two  shock  waves  (traveling 
along  x/at  =  ±1)  and  a  contact  discontinuity  at  a:  =  0.  The  bar  remains 
in  the  low  strain  phase  at  all  times.  One  seeks  a  self-similar  solution  of  the 
form 

for  — oo  <  X  <  —at  , 

for  —at  <  X  <0  , 

ior  0  <  X  <  at  ,  ,  ^ 

(42) 

for  at  <  X  <  oo  , 

where  are  the  only  unknowns.  Notice  that  the  continuity  of  T  across 

the  shock  waves  and  the  continuity  of  6  and  v  across  the  contact  discontinuity 


S{x,t)  ,  v{x,t)  ,  T{x,t)  =  < 


Sl  iVL  ,Ti 

(5o  ,vo  ,Ti 
So  ,  Vo  ,  Tr 
Sr  ,vr  ,Tr 
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were  enforced.  The  solution  is  shown  in  Figure  3.  Let 


6r  +  6l  vr  —  vl 


and 


j  = 


VR  +  VL  . 
2 


then,  from  the  jump  condition  (38)  it  follows  that 


So  =  h  , 
Vo=j- 


(43) 


(44) 


It  was  assumed  that  the  bar  was  initially  in  the  low  strain  phase,  then,  for 
the  whole  bar  to  remain  in  the  low  strain  phase  at  all  times  and  from  (42),  it 
is  required  that  -1  <  (5o  <  min  [(5M(rL),JM(Tk)].  Hence,  from  (44),  there  is 
a  solution  of  the  form  (42)  to  the  Riemann  problem  if  and  only  if  the  initial 
datum  h  is  such  that 


— 1  <  h  <  min  [5m{Ti,),Sm{Tr)]  .  (45) 


5.2  Solutions  with  two  phase  boundaries. 

The  solution  with  two  phase  boundaries  involves  two  shock  waves  (traveling 
along  x/at  =  ±1),  a.  contact  discontinuity  at  x  =  0  and  two  phase  boundaries 
at  x/ at  =  v  <  0  and  x/at  =  v,  >  0.  One  seeks  a  self-similar  solution  of  the 
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form 


S(x,t)  ,  v(x,t)  ,  T{x,t) 


f 


< 


V 


4  ,  ,  Ti  for 

5o  ,vo  ,Tl  for 

(5i  ,vi  ,T'  for 

(5i  ,iii  ,T"  for 

^2  ,V2  ,  Tr  for 

Sr  ,vr  ,  Tr  for 


— oo  <  X  <  —at 
—at  <  X  <  awt  , 
avt  <  X  <  0  , 

0  <  X  <  av,t  , 
av,t  <  X  <  at  , 
at  <  X  <  oo  , 


(46) 


where  Sq,  S2,  vo,  Vi,  V2,  T'  and  T"  are  unknown.  Notice  that  the  continuity 
of  T  across  the  shock  waves  and  the  continuity  of  5  and  v  across  the  contact 
discontinuity  were  enforced.  Here,  as  in  the  previous  section,  the  bar  is 
initially  in  the  low  strain  phase.  The  form  of  the  solution  is  shown  in  Figure  4. 
Prom  the  jump  conditions  (38),  (39)i,2  and  (40)i,2,  one  has 


-(<5o  -  Sl)  +  (uo  -  ul)  =  0  , 

v(5i  —  ^o)  +  (^1  ~  t'o)  =  0  , 

v,((52  -  ^1)  +  {V2  -  Ui)  =  0  , 

(<^i?  ~  S2)  +  (vr  —  V2)  =  0  , 

(^1  -  Jo  -  1)  +  (Ui  -  t^o)v  =  0  , 
(S2  —  Ji  +  1)  +  (V2  —  t5i)v,  =  0  . 
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Solving  this  system  provides 


and 


<^1  —  h  -\- 


1 


Vo 

Vi 

V2 


1 

1 

,1  +  v. 

1  +  v 

'  1 

1 

1 

1  —  V 

1  +  v* 

'  1 

1 

1  -V 

1-V* 

1 

1  ’ 

1  +  v* 

1  +  v 

1 

1  ' 

1  +  V* 

1-v 

1 

1 

1-V. 

1  —  V 

(47) 


(48) 


Moreover,  using  the  jump  conditions  (39)3  and  (40)$  together  with  (46)  yields 


T'  —  -  (5i  +  So  —  1)  +  It  +  Tl  , 
T"  =  -  {S2  +  Si  —  1)  +  It  +  Tr  . 


Furthermore,  using  (47)  in  the  above  equations  provides 


r'  =  Ti  +  r(v,v,)  , 

(49) 

r"  =  rfl  +  r*(v,v*) , 

(50) 
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where 


Notice  that 

r*(-v*,-v)  =r(v,v*)  . 

It  is  important  to  remark  that  a  solution  of  the  form  (46)  is  not  unique.  In 
fact,  equations  (47)-(52)  represent  a  two-parameter  family  of  solutions  where 
the  phase  boundary  velocities  v  and  v*  can  be  used  as  such  parameters.  The 
uniqueness  issue  is  addressed  in  Section  6  with  the  postulation  of  a  kinetic 
relation  which  singles  out  an  appropriate  pair  (v,  v*)  to  fully  determine  the 
solution  with  two  phase  boundaries.  A  preliminary  analysis  regarding  the 
existence  of  a  solution  is  developed  in  the  next  section. 

5.3  Phase  segregation  conditions  and  entropy  inequal¬ 
ity. 

Prom  (45)  it  is  clear  that  for  suitable  initial  data  there  exists  a  solution  to 
the  Riemann  problem  with  no  phase  boundary.  The  solution  which  involves 
two  phase  boundaries  requires  further  analysis.  One  has  to  ensure  that  the 
solution  satisfies  the  entropy  inequalities  (39)4  and  (40)4  and  that  (5o,  Tl)  G 
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Pi,  iSi,r)  €  Pz,  {5i,T")  e  P3  and  {82, Tr)  €  Pi,  i.e. 


>  5m(r') , 

>  5m  (T")  , 

8q  e  (— 1,5m(Pl)]  , 
82  G  (— l,5M(T/i)]  , 

P'G(o,rc) , 

T"  e  {0,Tc)  . 


(53) 


Conditions  (53)  are  referred  to  as  the  phase  segregation  conditions.  Prom 
(47),  (49),  (50),  (51),  and  (28),  it  follows  that  the  restrictions  (53)  and  the 
entropy  inequalities  (39)4  and  (40)4  are  equivalent  to 


Po(Pl,v,v*)  <  (l-m)/i, 

Pi(Ti{,v,v,)  <  (l-m)/i, 

ff2(v,v.)  <  h  <  P3(Pl,v,v.)  , 

P4(v,v*)  <  h  <  H5{Tr,\/,\/,)  , 

P6(Pi,v,v.)  <  h  <  H7{Tl,V,V»)  , 

Hs{TR,\/,Vt)  <  h  <  Hg{TR,v,y,)  , 
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where  the  functions  Hi  are  defined  for  —  1  <  v  <  0,  0  <  v,  <  1  by 


Ho(ri,  V,  V.)  =  ic  +  m(Tt  -  To)  -  i(l  -  m)  I j  ^ 


m 


H2{v,\f*) 

H3{Tl,V,V,) 

i/4(v,V*) 

H5(Tr,\/,V,) 

v,v*) 

H7{Tl,V,V,) 

H8{TR,y,v^) 

Hg{TR,v,y^) 


1-] 

Fv*/ 

fc  +  m(r*-To)-i(l-m){Y^  +  Y:^} 

f  {l  +  Y^I+m/r. 

_ ^\-l 

2\l+v  1+vJ 

<5c  +  1  +  M(Tx,  —  Tc)  +  H2{y,  V*)  , 

1/^ _ L\_l 

2\l-v,  l-vj 

5c  +  1  +  M(Th  —  Tc)  +  H4{\i,  V,)  , 


21 


Notice  the  following  symmetries: 


■ff2(v,v*)  =  i74(-v.,-v)  , 

-f^3(-,v,v.)  =  Hsi-,  -V*,  -v)  , 

■f^6(-,v,v,)  =  Hsi-,  -V*,  -v)  , 

■f^7(-,V,V*)  =  i/9(-,  -v„  -v)  . 

Equations  (54)i,2,3,4  are  equivalent  to  (53)i,2,3,4  respectively;  the  upper  bounds 
in  (54)5,6  are  equivalent  to  the  upper  admissible  value  of  the  temperature  in 
(53)5,6,  whereas  the  lower  bounds  in  (54)5,6  are  equivalent  to  the  entropy  in¬ 
equalities  (39)4  and  (40)4  respectively.  Notice  that  if  the  entropy  inequalities 
are  satisfied,  then  T'  —  Ti  =  r  >  0  at  the  phase  boundary  x/at  =  v  and 
T"  —  Tjr  =  r*  >  0  at  the  phase  boundary  x/at  =  v*.  Thus,  necessarily,  the 
conditions  T'  >  0  and  T"  >  0  are  satisfied. 

The  inequalities  (54)  determine,  for  fixed  Tl,  Tr,  a  region  in  the  (v,  v*,  h)- 
space  where  it  is  possible  to  have  a  two-parameter  family  of  solutions  of  the 
form  (46).  For  simplicity,  consider  the  case  when  It  =  0  which  corresponds 
to  twinning.  Furthermore,  since  it  was  assumed  that  a  =  0,  then  (30)  and 
(33)  give 

M  =  -m  >  0  ,  ^c  =  ^,  Sc>  MTc  -  (55) 
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Under  these  assumptions,  the  functions  Hi  to  Hg  take  the  following  form: 


Wo(Ti,v,v.)  =  M(Tc  -Ti)  -  1(1  +  M)  I  -  l| 


M  r  1 

'''  T  1  1  -  V2 


/f.(TR,v,v.)  =  M(rc  -  Ts)  -  1(1  +  M)  I  -  ll 

M  r  1  1 


1  +  V  1  +  V, 


Jf3(Ti,v, V.)  =  --M(Tc-Ti)  +  i/2(v,v.)  , 

H^(Tr,  V,  V.)  =  5  -  M(Tc  -  Tr)  +  H,(M,y.), 
«.(v,v.)  =  i/2{v,v.)  +  ?-i{j^)  , 
HriTL,^,^*)  =  {Tc  —  Ti}  +  if6(v,v»)  , 


■f^8(v,v*)  =  i/4(v,v.)  + 


3  1  f  1 


2  2  1  -  V 


HgiTfi,  V,  V*)  =  {Tc  —  Tr}  +  Hs{v,  v,)  . 


TT  _  '^0  tT  _ 

“  1  +  M  ’  ^  “  1  +  M 


Notice  that 


m  >  Hg  if  v"  > 


2  ^  (MTc  -  1/2)  -  MTl 


M{Tc-Tl) 

But,  from  (55)3,  {MTc  - 1/2)  -  MTr  <  0,  hence  Hg  >  Hg  for  all  v.  A  similar 


analysis  reveals  that  >  Hi  for  all  v*.  Moreover,  one  has 


H2 

> 

Ha 

if 

< 

V^ 

He 

> 

He 

if 

< 

V^ 

H2 

> 

He 

if 

v2 

> 

Cl  , 

H4 

> 

He 

if 

V2 

* 

> 

cf, 

where 


The  lower  bound  for  h  is  given  by 


for  v2  >  Hmin  = 


He  for  0  <  <  cf  , 

H2  for  cf  <  v^  <  1  , 


for  vl  >  v2  Hrt,in  = 


He  for  0  <  V*  <  cj  , 
H4  for  <  1  . 


Henceforth,  for  definiteness,  it  is  assumed  that 


(56) 


TR>n.  (57) 

For  the  upper  bound,  one  has  to  consider  two  cases  for  the  given  initial 
temperatures  Tl  and  Tr.  Only  the  first  case  is  considered  here:  Assume  the 
given  temperatures  Tr  and  Tr  satisfy 

2(M  +  1)(Th-Tx,)<1.  (58) 
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Under  this  assumption,  it  follows  that  Hy  <  Hz  and  <  Hz  for  all  and 
v^.  Moreover, 


Hmax  — 


CD 

for 

v^  <  1  - 

Hy 

for 

>  1- 

1 

4(TH-ri,)  +  l/(l-v2)  ’ 
4{Tr-Tl)  +  1/{1-w^)  • 


(59) 


The  lower  and  upper  limits  for  h  in  the  v^,v^-plane  corresponding  to  the 
phase  segregation  conditions  are  shown  in  Figure  5.  Notice  that  Hg  >  Hg 
for  any  v^,v*,  hence  the  admissible  region  is  always  non-empty.  The  values 
of  v^,  V*  for  which  iifmax  >  Hmin  Correspond  to  the  shaded  region  in  Figure  6. 
Thus,  the  shaded  region  corresponds  to  the  projection  on  the  v^,v*-plajie 
of  the  admissible  region  in  the  v^,  v^,  /i-space.  Typical  cross  sections  of  this 
region  in  the  v^,  h  and  v*,  h-planes  are  shown  at  the  end  of  next  section  in 
Figures  7  and  8. 


6  Driving  traction,  kinetic  relation  and  nu- 
cleation  criterion. 

As  pointed  out  in  Section  5.2,  the  solution  involving  two  phase  boundaries 
is  not  unique.  The  postulation  of  a  kinetic  relation  at  each  phase  boundary 
settles  the  uniqueness  issue.  In  order  to  introduce  the  kinetic  relation,  the 
explicit  form  of  the  driving  traction  acting  on  eaeh  phase  boundary  is  given 
below. 
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6.1  Driving  traction  for  phase  boundaries. 

For  the  thermoelastic  material  undergoing  an  adiabatic  process,  equation 
(26)  at  the  phase  boundary  which  has  the  high  strain  phase  on  its  right 
becomes,  with  a  =  0, 

.  (60) 

If  the  high  strain  phase  is  on  the  left  of  the  phase  boundary,  then  (26) 
becomes 


/  =  pow{log^-^}  . 

Consider  a  dimensionless  driving  traction  defined  as 


P7t  ' 


(61) 


(62) 


From  (33),  (60),  (61),  (49)  and  (50),  the  driving  traction  at  the  phase  bound¬ 
ary  x/at  =  V  is 


f=-f|ri,r(»,v.)]  , 

and  the  driving  traction  at  the  phase  boundary  x/at  =  v»  is 

f.  =  F|rs,r.(v,v.)l  , 


(63) 


(64) 
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where 


i^(r.r)=(r+0{log(l  +  ^)-|}  .  (66) 

and  r,  r*  are  given  by  (51)-(52).  In  the  purely  mechanical  case,  Abeyaratne 
&  Knowles  [2]  showed  that,  under  certain  assumptions,  the  corresponding 
Riemann  problem  was  symmetric  in  the  sense  that  v  =  — v*  and  f  =  — f* 
so  that  the  problem  reduces  to  a  one-parameter  family  of  solutions.  The 
adiabatic  Riemann  problem  does  not  exhibit  such  symmetry.  One  has  to 
consider,  separately,  a  kinetic  relation  for  each  propagating  phase  boundary. 
Nevertheless,  one  can  easily  check  that  the  solution  is  symmetric  when  Tl  = 

Tr. 

6.2  Kinetic  relations. 

From  (19),  (21),  (22)  and  (23),  the  dimensionless  rate  of  entropy  production 
for  a  segment  of  the  bar  containing  a  phase  boundary  during  some  time 
interval  can  be  expressed  by 

=  1  fv 

r  =  — r  =  —  . 

pac  (T) 

If,  following  [1]  and  [3],  the  quantity  f/{T)  is  identified  as  an  affinity  and  v 
as  the  corresponding  flvx,  then  the  kinetic  relation  relates  the  affinity  to  the 
flux  and,  in  this  case,  to  the  temperature  on  one  side  of  the  phase  boundary. 
The  temperature  is  not  considered  as  an  affinity  (it  describes  the  state  of  the 
material).  Therefore,  the  relation  between  the  affinity  and  the  flux  is  of  the 
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form 


Tl  +  ’'/2 
f. 


TR  +  rj2 


=  (/?13(TL)V)  , 


Clearly,  it  is  required  that 

/  ^  (69) 

</^13(-,v)  =  V). 

The  hine«c  relation  has  to  satisfy  sonre  restrictions 

olUv  Since  V  is  non-negative,  then  (/?ij 
ment  of  the  entropy  inequality,  bin 

such  that 

<Pij(r,v)v>o. 

This  implies  that  if  is  differentiable  with  respect  to  v.  then 

..i(r.o)  =  o.  ?^(T.0)>0.  P®) 

.n  .  is  strictly  monotonically  increasing  with  v. 
Furthermore,  assume  t  a  ¥>„  (67),  (68),  (69),  (63), 

Returning  to  the  special  case  where  h  -  0.  from 

(64)  and  (65),  one  has 

,og(l  +  ^)=^3r(Tt.-v).  log(l  +  ^)=^s.(Ts,v.)i 
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therefore, 

r  =  ,  r,  =  -  l)  . 

Using  (51)  and  (52)  with  It  =  0  provides 

(71) 

(72) 

Consider  the  restriction  of  the  function  <pz\{T,  •)  to  non-negative  values  of  its 
argument.  Define 

^31  {T,  v^)  =  ipzi  (T,  v)  for  V  >  0  .  (73) 

Recall  that  the  phase  velocities  axe  — v  >  0  and  v*  >  0,  hence,  subtracting 
(72)  from  (71)  and  using  (73)  gives 

+  2Tr  -  l)  ■ 

(74) 

Prom  the  properties  of  the  kinetic  relation  (pzii  the  function  (pzi  increases 
monotonically  with  v^.  As  a  function  of  v^,  (resp.  v*),  the  left-hand  side 
(right-hand  side)  of  (74)  is  strictly  monotonically  increasing  from  0  at  =  0 
(v^  =  0)  to  -t-oo  at  =  1  (v^  =  1).  Thus,  there  is  a  unique  €  [0, 1)  for 
a  given  e  [0,1).  Therefore,  it  is  possible  to  define  a  functional  relation 
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between  and  v^,  say 


v^  =  ^^,lTL,TR)  .  (75) 

Notice  that,  since  v  <  0,  then  v  =  - .  Henceforth,  let 

=  <f3l{TR,vl)  . 


Adding  (71)  and  (72)  and  using  (73),  (75)  and  the  above  notation,  one  has 


h  =  'i{vlTL,TR) 

/  V/  1  \  ^  WB  1  \  I  ^  f 


2  ^  2V-V^ 


(76) 


The  function  ^  is  strictly  monotonically  increasing  in  from  0  at  v*  =  0  to 
+00  at  =  1.  Therefore,  for  a  given  h  in  the  admissible  range,  (76)  singles 
out  a  unique  which,  by  (75),  provides  a  unique  v^.  This  fully  determines 
a  unique  solution  of  the  form  (46). 

The  fact  that  a  solution  exists  can  be  seen  by  the  following  analysis.  The 
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lower  bounds  Hq  and  Hs  of  the  admissible  region  can  be  expressed  as 


H6{<,Tl,Tr)  =  Hemvln,TR),yl) 

_1  fl  +  V-V^  1  ) 

~  2  t  1-V^  1  +  v.j  ’ 

H8{<,Tl,Tr)  =  HsmsflTL,TR),wl) 

_  1  fl+V.-v2  1  ] 

“2\  l-v2  1  +  Vj  ' 

Notice  that  >  max^2[He{\/l,Ti,TR),Hsi'^l,TL,TR)].  Moreover, 

the  upper  bound  Hg  is  such  that  Hg  =  {Tc  —  Tr)  for  v*  =  v  =  0.  Hence, 
by  continuity,  there  exists  a  range  of  values  of  h  in  the  admissible  region  for 
which  max[i?65  Hs]  <  h  =  '^  <  Hg,  thus  the  existence  of  a  solution  of  the 
form  (46)  is  guaranteed.  A  complete  analysis  of  solvability  requires  further 
knowledge  of  a  specific  kinetic  relation. 

6.3  Nucleation  criterion. 

As  shown  in  Sections  5.1  and  5.2,  there  are  two  different  kinds  of  solutions  for 
the  adiabatic  Riemann  problem.  One  can  interpret  the  no-phase  boundary 
solution  as  the  limit  of  the  two-phase  boundary  solution  when  v  ^  0~  and 
V*  —¥  O"''.  Now,  if  the  initial  datum  satisfies  (45),  it  is  possible  to  have  a 
solution  with  no  phase  boundaries  of  the  form  (42).  Assuming  that  Tr  >  Tr 
and  using  (55)i,  then  (44)  becomes 

-l<h<^-M{Tc-TR)  .  (77) 
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On  the  other  hand,  the  lower  bound  for  h  corresponding  to  the  solution  with 
two  phase  boundaries  at  =  0  is  /fmin(0, 0)  =  0  and  the  upper  bound  is 

■^^max(O)O)  =  {Tc  —  Tr).  Moreover,  the  curve  defined  by  (76)  passes  through 
the  origin,  hence,  if 


i  -  M(rc  -  Ti)  <  {Tc  -  Tfi) 

then  there  exists  an  overlapping  range  for  the  initial  datum  h  for  which  it  is 
possible  to  have  either  kind  of  solution: 

{-1  <  h  <  0 

,  then 

0  <  h  <  -  -  M(Tc  -  n) 

To  select  the  appropriate  solution,  a  nucleation  criterion  is  required.  Fol¬ 
lowing  [2],  assume  there  is  a  critical  value  fcr(T)  of  the  driving  traction  at 
which  a  transformation  from  low  to  high  strain  phase  occurs.  In  this  case, 
the  critical  value  depends  on  the  temperature  of  the  corresponding  low  strain 
phase  temperature.  From  (51),  (52),  (63),  (64)  and  (65),  nucleation  occurs 
when 


No  phase  change  solution. 
Both  types  of  solutions. 


(78) 


|f  I  =  F{Tc,  r(/i,  V,  V,))  =  f„(Ti)  (79) 


for  the  leftward  moving  phase  boundary  and 


f,  =  F(rH,r.(h,v,v*))  =  f„(rH)  (80) 
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for  the  rightward  moving  phase  boundary.  Although  h  is  considered  as  fixed 
(initial  conditions),  the  explicit  dependence  of  r  and  r»  on  h  is  shown  for 
clarity.  Equations  (79)  and  (80)  define  two  surfaces  on  the  v^,  v*,  h-space. 
The  intersection  of  the  region  where  f  >  fcr(7L)  and  |f*|  >  fcr(TR)  corresponds 
to  the  region  where  nucleation  occurs.  Typical  cross  sections  of  this  region 
are  shown  by  the  shaded  area  in  Figures  7  and  8. 

7  Nucleation:  nontrivial  solution  for  uniform 
data. 

The  importance  of  the  Riemann  problem  with  uniform  initial  data  is  related 
to  its  connection  with  the  nucleation  criterion.  Since  the  bar  is  initially  in 
the  same  phase  and  there  are  no  discontinuities,  the  existence  of  a  nontriv¬ 
ial  solution  (where  two  phase  boundaries  nucleate  from  an  arbitrary  point) 
could,  in  principle,  provide  a  restriction  on  the  critical  value  of  the  driving 
traction  f^.  Suppose  that  the  initial  data  of  the  Riemann  problem  of  Sec¬ 
tion  5.2  is  such  that  5l  =  5r  =  and  Ti  =  Tr  =  T^.  Since  the  system 
of  equations  (34)-(36)  and  the  jump  conditions  (38)-(40)  are  homogeneous, 
then  the  problem  admits  a  trivial  solution  S{x,  t)  =  6,  and  T(x,  t)  =  T»  Vx,  t. 
Nevertheless,  there  is  also  a  nontrivial  solution.  Since  Tr  =  Tr  and  in  view 
of  the  properties  of  the  kinetic  relation  (pij,  equation  (74)  implies  that 

=  v^  =>  V  =  — v*  . 
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Prom  (47)-(52),  there  exists  a  one-parameter  family  of  nontrivial  solutions 
for  the  Riemann  problem  with  uniform  initial  data  of  the  following  form 
(parametrized  by  v»): 


and 


where 


S{x,t) 


5*  —00  <x  <  —at  , 

So  —at  <  X  <  —aWtt  , 
*  5i  —av,t  <  X  <  avj,  , 
^0  <  X  <  at  , 

5,  at  <  X  <  oo  , 


r* 


T{x,t)  =  { 


r 


r. 


— oo  <  X  <  —aVtt  , 
— av*t  <x  <  omJ,  , 
av*t  <  X  <  oo  , 


(81) 


(82) 


Notice  that,  in  view  of  the  analysis  of  Section  6.2,  a  solution  of  the  form  (81), 
(82)  exists.  Using  the  notation  of  Section  6.2,  it  follows  that  U  =  v*  and 


34 


‘Pl  =  ^R  =  'P*  =  p{T*,yD-  Equation  (76),  together  with  (43),  becomes 


The  nucleation  criterion  was  introduced  in  terms  of  the  driving  traction 
which,  in  this  case,  can  be  written  as 


f(T.,r((5*,v.)) 


where 

In  order  to  have  an  admissible  solution,  the  phase  segregation  conditions  and 
entropy  inequality  requirements  have  to  be  enforced.  From  Section  5.3,  the 
lower  bound  for  the  value  of  h  is  given  by 


2v«  -  vj 

1  —  v2  _ 


Hmin  — 


0  <  <  cf  , 

cf  <  <  1  . 


Recall  that  the  functions  and  Hs  are  related  to  the  entropy  inequality 
(39)4  and  (40)4,  whereas  the  functions  H2  and  H4  represent  the  restriction 
T'  =  T"  <  Tc-  The  upper  bound  for  the  admissible  region  is  given  by 


The  admissible  region  is  shown  in  Figure  9.  FVom  this,  one  can  see  that  there 
is  a  minimum  value  for  below  which  no  nucleation  occurs  (in  this  case. 
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=  0  at  V*  =  0).  Observe  that  this  lower  bound  is  given  by  the  entropy 
inequality  (it  corresponds  to  f  =  0).  This  analysis  confirms  the  fact  that 
nucleation  occurs  at  a  critical  value  such  that  fcrV*  >  0  (as  required  by  the 
entropy  inequality)  but  fails  to  restrict  the  nucleation  criterion  in  any  other 
way. 
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Low  strain  phase  Unstable  phase  High  strain  phase 


Figure  1;  Material  phases  in  the  temperature-strain  plane. 
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Figure  3:  Solution  with  no  phase  boundary. 
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Figure  4:  Solution  with  two  phase  boundaries. 
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